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This is a brief progress report on a long-term collaborative project to measure the power 
spectrum (PS) of nws density fluctuations from the Mark III and the SFI catalogs of 
peculiar velocities.Ba The PS is estimated by applying maximum likelihood analysis, 
using generalized CDM models with and without COBE normalization. The applica- 
tion to both catalogs yields fairly similar results for the PS, and the robust results are 
presented. 

1 Introduction 

In the standard picture of cosmology, structure evolved from small density fluctua- 
tions that grew by gravitational instability. These initial fluctuations are assumed 
to have a Gaussian distribution characterized by the PS. On large scales, the fluc- 
tuations are linear even at late times and still governed by the initial PS. The PS 
is thus a useful statistic for large-scale structure, providing constraints on cosmol- 
ogy and theories of structure forniatioiL In recent years, the galaxy PS has been 
estimated frorru^everal redshift surveys.u In this work, we develop and apply like- 
lihood analysigj in order to estimate the mass PS from peculiar velocity catalogs. 
Two such catalogs are used. One is the Mark III catalog of peculiar velocities,El 
a compilation of several data sets, consisting of roughly 3000 spiral and elliptical 
galaxies within a volume of ~ 80h~^Mpc around the local grjoup, grouped into 
~ 1200 objects. The other is the recently completed SFI catalog,Q a homogeneously 
selected sample of ~ 1300 spiral field galaxies, which complies with well-defined 
criteria. It is interesting to compare the results of the two catalogs, escecially in 
view of apparent discrepancies in the appearance of the velocity fields.QEI 

2 Method 



Given a data set d, the goal is to estimate the most likely model m. Invoking 
a Bayesian approach, this can be turned to maximizing the likelihood function 
£ = ^(dlm), the probability of the data given the model, as a function of the 
model parameters. Under the assumption that both the underlying velocities and 
the observational errors are Gaussian random fields, the likelihood function can be 

written as £ = [{2tt)'^ det{R)]-^^^ exp (^-^J2^j diR'/dj^j , where {d,}f^^ is the 
set of observed peculiar velocities and R is their correlation matrix. R involves the 
theoretical correlation, calculated in linear theory for each assumed cosmological 
model, and the estimated covariance of the errors. 

The likelihood analysis is performed by choosing some parametric functional 
form for the PS. For each assumed PS, one can calculate the likelihood function 
and going over the parameter space find the the PS parameters that provide the 
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Figure 1: Likelihood analysis results for the flat ACDM model with h = 0.6. InC contours in 
the Q — n plane are shown for SFI (left panel) and Mark III (middle). The best-fit parameters 
are marked by 's' and 'm' on both, for SFI and Mark III respectively. The right panel shows the 
corresponding PS for the SFI case (solid line) and for Mark III (dashed). The shaded region is 
the SFI 90% confidence region. The three dots are the PS calculated from Mark III by Kolatt and 
Dekel (1997),^° together with their la error- bar. 

maximum likelihood. Confidence levels are estimated by approximating —2hiC as a 
distribution with respect to the model parameters. Note that this method, based 
on peculiar velocities, essentially measures /(r2)^P(fc) and not the mass density PS 
by itself. Careful testing of the method was done using realistic mock catalogs,cl 
designed to mimic in detail the real catalogs. 

We use several models for the PS. One of these is the so-called F model, where 
we vary the amplitude and the shape-parameter F. The main analysis is done with a 
suit of generalized CDM models, normalized by the COBE 4-yr data. These include 
open models, flat models with a cosmological constant and tilted models with or 
without a tensor component. The free parameters are then the density parameter 
ri, the Hubble parameter h and the power index n. The recovered PS is sensitive 
to the assumed observational errors, that go as well into R. We extend the method 
such that also the magnitude of these errors is determined by the likelihood analysis, 
by adding free parameters that govern a global change of the assumed errors, in 
addition to modeling the PS. We find, for both catalogs, a good agreement with the 
original error estimates, thus allowing for a more reliable recovery of the PS. 

3 Results 

Figure |l| shows, as a typical example, the results for the flat ACDM family of models, 
with a tensor component in the initial fluctuations, when setting h ~ 0.6 and varying 
ri and n. The left panel shows the \nC contours for the SFI catalog and the middle 
panel the results for Mark III. As can be seen from the elongated contours, what 
is determined well is not a specific point but a high likelihood ridge, constraining a 
degenerate combination of the parameters of the form rtn^ '^ = 0.59 ± 0.08, in this 
case. The right panel shows the corresponding maximum-likelihood PS for the two 
catalogs, where the shaded region represents the 90% confidence region obtained 
from the SFI high-likelihood ridge. 

These results are representative for all other PS models we tried. For each 
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catalog, the different models yield similar best-fit PS, falling well within each oth- 
ers formal uncertainties and agreeing especially well on intermediate scales (fc ~ 
0.1/iMpc"i). The similarity, seen in the figure, of the PS obtained from SFI to 
that of Mark III is illustrative for the other models as well. This indicates that 
the peculiar velocities measured by the two data sets, with their respective error 
estimates, are consistent with arising from the same underlying mass density PS. 
Note also the agreement with an independent measure of the PS from the Mark III 
catalog, using the smoothed density field recovered by POTENT (the three dots).E3 
The robust result, for both catalogs and all models, is a relatively high PS, with 
P{k)n^-^ = (4.5 ± 2.0) X 10^ (/i^iMpc)3 at A: = 0.1 /iMpc^i. An extrapolation to 
smaller scales using the different CDM models gives (Tgil^'® = 0.85 ±0.2. The error- 
bars are crude, reflecting the 90% formal likelihood uncertainty for each model, the 
variance among different models and between catalogs. The general constraint of 
the high likelihood ridges is of the sort fl h^o'^ n'^ = 0.75 ± 0.25, where fj, = 1.3 and 
v = 3.7, 2.0 for ACDM models with and without tensor fluctuations respectively. 
For open CDM, without tensor fluctuations, the powers are /i — 0.9 and ly — 1.4. For 
the span of models checked, the PS peak is in the range 0.02 < k < 0.06/iMpc~^. 
The shape parameter of the F model is only weakly constrained to F = 0.4 ± 0.2. 
We caution, however, that these results are as yet preliminary, and might depend 
on the accuracy of the error estimates and on the exact impact of non-linearities.a 
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